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We analyze the disturbed flow in the supersonic laminar boundary layer when local 
heated elements are placed on the surface. It is exhibited that these flows are described 
in terms of free interaction theory for specific sizes of thermal sources. We construct the 
numerical solution for flat supersonic problem in the viscous asymptotic layer in which 
the flow is decribed by nonlinear equations for vorticity, temperature with the interaction 
condition which provides influence of perturbations to the pressure in the main order. 



1. Introduction 

Th e present work continues the studies started in ( Lipatov 2006t Koroteev fc Lipatov 
20091) . which were devoted to construction of asymptotic solutions of Navier- Stokes equa- 
tions in the regions, in which local heating elements are situated on the surface of the 
body. In iKoroteev fc Lipatovl (|2009t ) we demonstrated how this problem can be solved 
analytically for small temperature perturbations. 

The idea lying behind this approach consists in the utilizing of small temperature 
perturbations for the control of separation of the boundary layer and delay of laminar- 
turbulent transition. The former question is related to the zero shear stress point in the 
flow which can alter its location if one affects the boundary layer by some, not necessarily 
thermal, perturbations. The latter is related to the possibility to slow down the flow by 
means of the same perturbations to decrease local Reinolds number Re. 

The main mechanism which enables to carry out the control of the boundaty layer 
is alteration of the pressure induced by that of displacement thickness in the flow from 
perturbations. The source of the perturbations, as we said, can have various nature, e.g., 
variation of curvature of the surface which can be realized as small humps or irregularities 
on it. The problems about perturbations in the boundary layer which are produced by 
small humps located on the surface in the botto m of the boundary layer served the object 
of thorough studies during p ast several decades ( Bogolepov fc Neiland 1971 : Smith 1974 
Bogolepov & Neiland''li7i) 

In (Lipatov 2006) the general description of problems emerging when local heating 
elements are located on the surface was given. In the same work an important similarity 
between the problems in question and those of flows over the surface with small humps 
was demonstrated. It was, in part, shown that the local heating forms an effective hump 
and the outline of the flow becomes sim ilar to that of physical humps on the surface 
(|Smithlll97j: iBogolepov fc Neilandlll976l ). 

The purpose of the present paper is to extend the analysis to more general and more 
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complicated case of nonlinear perturbations. Linear problems, corresponding to small 
perturbations AT << 1, though being novel, only presented the first step in studing flows 
in the boundary layer with thermal pertu rbations. In addition, the nonlin ear problem 
studied here, unlike the linear formulations ( Koroteev fc Lipatovll2009l 2011 ). requires an 
extensive use of numerical methods and usually nontrivial algorithms for computation. 



2. Formulation of problem for the viscous sublayer 

We consider the uniform supersonic flow over a fiat semi-infinite plate for large Reinolds 



numbers Re when laminar-turbulent transition is absent. Denote Re 



and consider 



it as a small parameter. The fiow is assumed to be fiat and stationary. On the surface of 
the body a heating element of size a is located which produces thermal perturbations in 
the fiow. The size of the element is an important parameter of the problem and discussed 
below. We assume that thermal perturbations are not small AT T ^ 0(1). 

We are interested in constructing asymptotic solutions of Navier-Stokes equations sup- 
plied in this case by the equation for temperature, in regions where local perturbations 
of the temperature exist. This construction is accomplished in terms of the interaction 
theory by divi ding the region, in the vicinity of the leating element by th ree smaller re- 
gions or decks (Neiland 19691 : Stewartson fc Williams 1969t Messiteij[l970[ ). Sizes of these 
regions are proportional to powers of the parameter Re = e~^. The role of these regions 
consists in different character of perturbations in them and different infiuence on the fiow 
as a whole. For the problem under cons ideration the analy sis of equations corresponding 
to each region was fulfilled in .Koroteev &: Lipatov ( 2009.) ; iLipatov ( 2006 ) and therefore 
is not described here. N ote that t his analysis shows th at the flow is described by the 
regime of free interaction iNeiland ( 1969[ ): Lipatovl ([2006'). In this regime thermal distur- 
bances induce pressure gradient which influences the fiow functions in the main order. 
That means that diffusive, convective, and press ure gradient te rms in the Navier-Stokes 
equations have the same order. As exhibited in ( Lipatov 2006 h this condition furnishes 
the following asymptotic relations Ap ^ Ofe^/^), o ^ 0(e^/^), Ay ~ 0(e^^'^) which cor- 
respo nd to the known scales of free interaction (Neiland 1969t Stewartson fc Williams! 
1969( ). Thus, the free interaction regime emerges only for specific sizes of the thermal 



region. 

From this analysis it also follows that the upper deck is described by inviscid Euler 
equations, the middle deck by Prandtl boundary layer equations and the lower layer, 
the closest one to the surface, by some nonlinear set of equations of parabolic type 
( Koroteev fc Lipatov 2009i I2OIII) . The asymptotic layout of three decks is portreyed in 
fig. [1] In addition, the equations on the lower deck are supplied by boundary and initial 
conditions, which are formulated below, and also by an interaction condition which not 
only connects functions of lower deck to those of upper deck, but also provides infiuence 
of the right boundary of the whole fiow region to the behavior of functions inside the 
region. It implies that the set of equations in the lower deck can not be viewed as that 
of parabolic type and gains, in a sense, the property of ellipticity. 

The lower deck is the most essential part of the hierarchy for the pressure gradient 
presented in equations for all decks is not fixed but becomes self- induced, i.e., varies on 
account of displacement thickness of the boundary layer to the outer inviscid flow, and 
the main contribution to the variation of the pressure is furnished by the very lower deck, 
while the main deck(the main part of the boundary layer) remains passive and merely 
conveys perturbations from the lower deck to the outer flow. The indicated scheme enables 
to construct solutions of equations in the lower deck with self-induced pressure gradient 
and thus provides interaction between lower and upper decks. 
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Figure 1. Sketch of triple deck region in the vicinity of heating element. 



We give the formulation of th e set of eq uations for lower deck, further called the viscous 
sublayer, referring to (Korotcc v fc Lipat ov 2009, 2011) for details of derivation. 
The set of equations for the viscous sublayer has the form 

dub ^ dvb_ 
dxb dyb 



0, 



dub dub 

"^b-^ ^ Vb-^ 

dxb oyb 



T,, 



dpb d^ub 



Ub 



dTb 



Vb- 



dxb 



(2.1) 



dxb ' ""dyb dyl 



Here Ub,Vb are longitudal and vertical components of velocity in the sublayer, Tb - the 
perturbation of temperature, pb{x) is the pressure which only depends on x and conse- 
quently is constant through the decks. 

The system is supplied by the following boundary conditions. On the surface the non- 
slip conditions are stated 

Ubixb,0) = Vbixb,0) ^0. (2.2) 
We prescribe a temperature profile on he surface: 

Tb{xb,0)^T^ixb). (2.3) 
The boundary condition at significant distance from the surface is presented as follows 



Ub 



r + txi 

yb + A{xb), Ai{xb) ^ I {I - Tb)dri + A{xb) 
Jo 



Tbix,(X}) ^ l,yb ^ (X). (2.4) 
The function A{xb) is a standard displacement by mean s of which the proble ms of in- 



terac tion without heat perturbations are formulated (Jobc fc Burgr afi 1974; Sv chev et 



1998). Local perturbations of the temperature generate an effective local hump ljLipatov 



20061) and additionally contribute to the boundary condition for j/f, — )■ -|-oo, which is 
expressed by means of a function Ai{xb)- The temprature perturbations vanish far from 
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the location of the heated region and A{xb) tends to Ai{xb) as it follows from the above 
condition. 

The conditions which relate the displacement thickness with the longitudal velocity are 
supplied by the interaction condition, which, in turn, relates the displacement thickness 
with the pressure which provides the interaction of the viscous sublayer with the outer 
inviscid flow. 

The condition of interaction differs for subsonic and supersonic flow s. We do not give 



its derivation for th is question was many times discussed in literature (|Stewartsonlll981 



Svchev et al. 1998() . In part, for the supersonic flow which we study, the condition of 



interaction has the form 

dAi{xb) , . 

Here i? is a constant which allows to consider different regimes of interaction. In our case 
B ^ 1 which corresponds to the regime of free interaction. 



3. Numerical solution of boundary problem 

The numerical solu tion the prob lem (|2.HI2.5p is constructed on the premise of the 
method suggested i n ( Ruban 19761) . the detailed description of which can be found in 



( Svchev et al. 19981 ). The method was originally employed for computations of subsonic 



problems but can be also utilized for supersonic flows. Supplying the equation for tem- 
perature we thus give but a short sketch of the procedure, paying attention to necessary 
modifications of the numerical method. 

In the first place the boundary problem described in the previous section is reduced 
to that for vorticity r = du/dyb- This is carried out by differentiating the momentum 
equation wrt. yb and taking into account the equation of continuity. Then the system 
becomes 



dr dr dTb dpb d^r 
dxb oyb ayb oxb oy^ 



dTb , dTb d^Tb 

Ub- — h Vb— — = ^r-r- 



(3.1) 



dxb dyb dy^ 
The boundary conditions have the following form 

T ^ l^yb +oo,Xb — s> — oo. (3-2) 



Non-slip conditions ()2.2|) as well as those for the temperature ()2.3|) are retained on the 
wall. 

Next it is easily seen that the momentum equation gives 

dp dT 

Tw{xb)^—^^, yb = 0- (3.3) 
dxb dyb 



From this equation we obtain using the interaction condition 

_ dT_ 
dxl dyb ' 

Localization of the effective hump implies the decay of perturbations downstream. 



BTn,{xb)^^ = —, 2/h = 0. (3.4) 



Local temperature perturbations in the boundary layer in regime of free viscous-inviscid interactioh 



namely 

p{xb) —J' 0, Xb ^ +00 
The last condition implies, according to (|2.5p . that 

dAi{xb) 



dxh 



-J- 0, Xb ^ +00. (3.5) 



The numerical method provides some freedom in the consecutive evaluations of the 
functions as the number of relations between them exceeds the number of functions. We 
describe but one possible realization of the procedure which enables to carry out the 
computations. 

The grid is given by {xj,yk), j — 1,2 . . . N, k ^ 1,2 . . . M, with steps Ax, Ay respec- 
tively. 

The system p.ip is approximated as follows 



dT\* I Tjk+l — Tjk-l I ^, Tjk+l—Tjk- 

jk 

_ Tjk-l — 2Tjk + Tjk+1 



2 Ay 2 Ay 



1 



Ay^ 



Tj+ik — Tjk , Tjk+1 — Tjk-i Tjk-1 — 2Tjk + Tjk+i 
A^ + 2A-y + A? 

where 

'dry J u,,^0 



(3.6) 



We thus notice that depending on the sign of the longitudal velocity different templates 
are emp loyed. This is quite a famous method for securing static stability of the difference 
scheme ( Roache 19721) . Der ivatives wrt. x can be approximated using the values in three 



points (|Svchev et al. [l99i). However our computations show that this alteration does 



not produce noticeable effect in accuracy, at least for supersonic flows. Therefore we 
employed tow-point templates. 

The solution is co nstructed by iterations with respect to Ai{xb) using under relaxation 
( Sychev et al. 1998h . Let the approximation of Ai{xb) is known on some step of the 



iteration procedure. At each line x — xj the first equation of (13. 6|) is reduced to a linear 
set of algebraic equations with a tridiagonal matrix and solved consecutively from the 
bottom to the top of the grid with the boundary conditions (|3.4p and (|3.2p and then 
the equation for temperature p.6p is solved with the boundary conditions (|2.3p . (|2.4p . 
Then the computations are transferred to the next line. It allows to fulfill computations 
of vorticity and temperature in the whole 2D region. 

The right boundary attained, the functions Ujk,Vjk as well as Ai{xb) are recomputed 
on the current iteration in the whole field and the new approximation for Ai (xb) is 
evaluated from the relaxation 

A'l+^ = {l-r)A'^+rAi{xb), 

where r is a relaxation parameter. In problems of interaction theory it usually tur ns out 
that convergence of itera tions is strongly sensitive to the values of this parameter! Smithl 



1977 : Svchev et al. 1998h . The problem under consideration is not an exception. In fig. [3] 



the errors AA = A^~^^ — A"' are presented for several values of the relaxation parameter. 
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Figure 2. Behavior or error A A in consecutive iterations for various values of the relaxation 

parameter r. 

pressure perturbation in the sublayer 




Figure 3. Pressure in the viscous sublayer obtained from the numerical solution of the 
problem. Small B decline pressure perturbations even if AT is not small. 



It is noticed that convergence, which exists for r — 0.02, vanishes even for r = 0.031. In 
our computations we took r = 0.02, the value which was found empirically and which 
yields an estimate of the upper boundary of some interval of convergence. 



4. Results of computations 

The key functions to evaluate are distributions of pressure and shear stress on the 
surface. As we said, the self-induced pressure is constant through the layers and thus its 
values remain valid in the inviscid flow for the fixed Xb- The shear stress in turn can be 
thought of as a general function characterizing forces in the flow. 

In fig. [3] it is exhibited the pressure in the viscous sublayer for heat humps, of a 
somewhat artificial form, which are given by 
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Figure 4. Pressure in the viscous sublayer obtained from the numerical solution of the problem 
compared to hnear solution which corr esponds to small temperature perturbations and was 
presented in (jKoroteev fc Lipatovl [20091 ') . For both cases the heated region is determined by 
l|4.1|l with the size 10. 



Here h the perturbation of the temperature and a the size of the heated region. For 
computations we take h — 0.3 — 0.4. 

This form of temperature distribution is simpler for analytic study and in the same 
time enables to represent adequately qualitative behavior of functions in the flow. The 
results are presented for various B. It is noticed that the pressure perturbations decrease 
for small B even if AT ~ This situation corresponds to compensated interaction 

regime which is discussed elsewhere. 

The comparison of computations for nonlinear and linear cases, the latter correspond- 
ing to small temperature perturbations, is represented in fig. |4l The difference is espe- 
cially noticeable in the transient region and downstream which are not quite adequately 
approximated by the linear theory. 

The flow moves from the left to the right. When approaching the heated region the 
positive pressue gradient produces aditional deceleration of the flow, thus delaying pos- 
sible laminar-turbulent transition. Note that there are two possible separation points 
in the flow, namely boundaries of the heated region upstream and downstream (fig. [5]). 
They have to be controlled if the separation is not desirable. In the transient region 
between these points and over the heating element the flow accelerates affected by the 
negative pressure gradient which simultaniously increases the shear stress. Finally per- 
turbations of the pressur decay Behind the heated region and the flow againg slows down. 
The parameter B from (|3.4|) influences the amplitude of pressure perturbations. As we 
mentioned before the regime of free interaction implies i? ~ 1 . 




(4.1) 
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Figure 5. Shear stress in the boundary layer. 



5. Conclusions 

The problem whose solution is presented in this paper presumabily can essentially con- 
tribute to possibility of control of boundary layer and general understading of influence 
of thermal perturbations to both subsonic and supersonic flows. Possible closest devel- 
opements of this work imply generalization of the obtained solutions to three dimensional 
case which is of course of more interest for applications. From the other hand, another 
direction would be the generalization of these results to nonstationary flow where it is 
essential to study nonlinear perturbations and their propagation. Both these problems 
are under consideration and will be published elsewhere. 
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